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Abstract Several ways are discussed how to control 
the Landau-Zener tunneling in the Wannier-Stark sys- 
tem. We focus on a realization of this system with inter- 
acting and noninteracting ultracold bosons. The tunnel- 
ing from the ground band to the continuum is shown 
to depend crucially on the initial condition and sys- 
tem parameters and, more interestingly, on added time- 
dependent disorder - noise - on the lattice beams. 
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1 Introduction 

Bloch oscillations, Landau-Zener (LZ) tunneling, and 
Wannicr Stark ladders [1,2,3, 4, 5,6,7,8,9, are fundamen- 
tal quantum effects occurring in a system of electrons 
moving in a periodic potential and subjected to a con- 
stant electric field. Due to complications such as impu- 
rities, lattice vibrations, and multiparticle interactions, 
clean observations of these effects have been difficult [TU] ■ 
In recent years, ultra-cold atoms and Bose-Einstein con- 
densates in optical lattices have been increasingly used 
to simulate solid state systems and the above mentioned 
phenomena p[7ll51[TTl[T2l [LllP11[L3 ] . Optical lattices are 
nowadays easy to realize in the laboratory, and their 
parameters can be perfectly controlled both statically 
and dynamically, which makes them attractive model 
systems for crystal lattices. More complicated poten- 
tials can be realized by adding further lattice beams [TfJI 
[T71[TMT§ll20ll21"] . In fact, by superimposing laser beams 
from different directions and with slightly different wave- 
lengths, it is possible to generate many different three- 
dimensional lattice geometries [T4HT5] . The question arises 



of how to control the dynamics of particles by quasiperi- 
odic potentials (possibly time-dependent or even stochas- 
tic ones). 

In this paper, we present results on the Wannicr 
Stark system realized with ultracold atoms, forming a 
Bose-Einstein condensate, in an optical lattice [6j[7l[8j 
IT2l[22ll23ll24ll251l26j . We compute the time dependence 
of the tunneling probability of the Bose-Einstein con- 
densate atoms out of the ground band in which they 
were originally prepared. By changing the initial condi- 
tion and the system parameters and introducing atom- 
atom interactions into the system, we are able to control 
the tunneling rate of the Bose-Einstein condensate to 
higher bands. Finally a controlled noise added to the 
system will be shown to be a further handle to engineer 
the interband tunneling. 



2 Landau Zener tunneling in optical lattices 

We study the temporal evolution of ultracold atoms 
loaded into a quasi one-dimensional optical lattice 
which can be a spatially periodic potential or a time- 
dependent stochastic potential, subjected to an addi- 
tional static force, in the presence of weak atom-atom 
interactions [55]. We use the following general form of 
three-dimensional Gross-Pitaevskii equation to model 
the temporal evolution of the atoms 



H = + W (r, t) + g 3D |#(r, t)\ 2 4 

W (r, t) = T4 rap (r) + V x (x) + V 2 (x, t) , 



Fx: 



(1) 



with lattice potentials along the longitudinal direction, 
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VAx) = aVs\rL Z — 



V 2 (x, t) = aV sin^ — 



(2a) 



(2b) 



2 

M is the mass of condensate atoms and F is the Stark 
force. V\(x) and V2(x,t) are spatially periodic potentials 
with incommensurate lattice spacings g?l and d L , respec- 
tively. As will be shown below around Eq. ^ , the noise 
has a tendency to average over the second lattice, and 
therefore the amplitudes of the two lattices should be 
comparable. For convenience we chose equal amplitudes 
aV and — 426 nm and d L = cZl(\/5 — l)/2 for the 
lattice constants. V2(x,t) is a time-dependent stochastic 
potential with a time-dependent stochastic phase <j>(t), 
which we will characterize further down in section 12.21 
The renormalization factor a is introduced to be able 
to compare the dynamics in the presence of the poten- 
tial given by Eqs. (|2a|l and (|2b[) to the dynamics of the 
"reference system" , i.e. the dynamics in the potential 
W (r) = Vtrap (r) + V sin 2 (^^j. a will be chosen in such 
a way that the following standard deviations are equal: 

(Hf)-Mf)>Jl 

= ((V cS (x)-(V cS (x)) x f) x , (3) 

where the effective potential V c s(x) will be defined in 
section 12.21 The average (•) is an integral over space 

for a sufficiently large L, i.e., (•) = 1/L J_ L , 2 -dx. The 
third term in the Hamiltonian is the nonlinearity, which 
makes the equation different from the Schrddinger equa- 
tion, ^(r, t) represents the condensate wave function and 
|!f(r,i)| 2 the local atomic density. g 3D = 47r? ^° sjVo is 
the coupling constant which is proportional to the scat- 
tering length a s and determines the strength of atom- 
atom interactions, where 4 ^ E a ^ ~ 2.45 x 10~ 21 m 3 , 
with a s = 53 x 10" 10 m and M = 1.44 x 10" 25 kg 
for rubidium 87. -/Vo is the number of atoms in the con- 
densate. The recoil energy Er — p 2 l /2M and the recoil 
momentum j»r = nh/d^ are the characteristic energy 
and momentum scales for our system. Moreover, we de- 
fine Vq — V/En and Fq — Fd^/En as dimensionless 
quantities in this energy unit. Experimentally, the ini- 
tial state is the relaxed condensate wave function pre- 
pared in the confining potential given by a harmonic trap 
Vtrap = \m (u}pP 2 + oj 2 x 2 ), with u; x <C w p for a quasi ID 
situation, and then loaded adiabatically into the optical 
lattice when the Stark force F equals zero. Then w x is 
either switched off or relaxed to a small value uj XtVe i and 
the Stark force F is simultaneously switched on to in- 
duce the dynamics. In the following sections we study 
the system in the presence/ absence of different terms of 
the above Hamiltonian. 



2. 1 Noise free case 
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Fig. 1 Time resolved survival probability of the BEC in the 
ground band, for a lattice depth Vo = 4 E-r, and a narrow 
width of the initial momentum distribution of the BEC cloud, 
Z\p/2pR, = 0.05 (trap frequency tj K = 2n x 50 Hz), and var- 
ious Stark forces. For non-RET condition: Fq = 1.07 (solid 
line) and Fo = 1.63 (dash-dotted line); for RET condition: 
Fo = 1.48 (dashed line). Inset: height of one step (of data as 
shown in the main panel) for fixed Vq = 4 and various 7. 
The step height as predicted by Eq. (T3| (solid line) in compar- 
ison with the numerically obtained values (filled circles). A 
significant deviation from the LZ prediction can be observed 
at values of — corresponding to RET conditions and marked 
by the vertical dashed lines. 



2.1.1 Linear case - Wannier-Stark problem For zero 
or small w Xjrc i (c.f., e.g., [28 ) and negligible atom- 
atom interactions, when Vzix, t) = and a = 1, the 
above Hamiltonian in Eq. (JTJ) describes the dynamics of 
atoms in a tilted periodic potential which is the well- 
known single-particle Wannier-Stark problem. Without 
the nonlinearity term, we can use the ID version of 
Eq. ([T]) for our simulations. In the presence of F, the 
quasimomentum of a condensate (initially prepared at 
the center of the Brillouin zone in the ground band) 
scans the lower band in an oscillating motion periodi- 
cally - so-called Bloch oscillations - with the Bloch pe- 
riod T B = 2h(Fd L )- 1 . At the edge of the Brillouin zone, 
where the gap between the ground and the first excited 
band AE of the F — system (increasing with Vo [2T]). 
acquires its minimum value, a tunneling of the conden- 
sate to the first excited energy band may occur. The 
tunnelled atoms escape from the system through succes- 
sive tunneling events across the much smaller band gaps 
between the upper bands. This phenomenon is known as 
the LZ tunneling. LZ theory predicts a decay rate 

PLZ = e~5, (4) 

where 7 is the adiabaticity parameter and 7 w [57] . 
In order to study the LZ prediction for our system we 
need to access the decay rate of the population from 
the ground band. In that respect, we compute the time- 
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dependence of the probability of the condensate to re- 
main in the ground band in which it has been initially 
prepared. Such a survival probability is best measured in 
momentum space. From the time-dependent momentum 
distribution we can determine P sur (t) by projection of 
the evolved state ^(p, t) on to the support of the initial 
state [28U30] 

/ dp y d Pz dp x \<F(p,t)\ 2 , (5) 

- OO «/ — OO J —Pc 

where p c is an ad hoc cut-off momentum. In our cal- 
culation we chose p c — 3p roc . Then Eq. ([S]) starts to 
measure the wave packet leaving the ground-band with 
an artificial delay of one Bloch period, since by the accel- 
eration theorem the average momentum is proportional 
to time (more details can be found in [28ll30] h A very 
nice step structure - local deviation from the overall ex- 
ponential decay - can be seen in the survival probabil- 
ity (see Figures Q] and [2]). Such step structures reflect 
the above mentioned phenomena, i.e., Bloch oscillations 
and LZ tunneling. One can see that the tunneling events 
from the ground band to the next band occur after each 
Bloch period when the wave packet is at the edge of the 
Brillouin zone. It is possible to control the LZ tunneling 
of the Bose-Einstein condensate from the ground band 
in the linear system and in the absence of atom- atom 
interactions by 

— changing the system parameters, such as external 
force -F , and the amplitude of the optical lattice Vq 
(see Figured]), 

— exploiting resonantly enhanced tunneling (RET) be- 
tween degenerate Wannier-Stark states at AE rj 
n x Fd^, with n being an integer number (see Fig- 
ure[T]) [9], 

— changing the initial condition by changing the trap 
frequencies and hence preparing the Bose-Einstein 
condensate with different widths Ap of its initial mo- 
mentum distribution (see Figure [2]) ■ 

As seen in Figure [TJ increasing the tilting force leads to 
more and more tunneling of the atoms from the ground 
band. On the other hand, depending on the system 
parameters, one can tune into a special condition for 
which the rate of the tunneling is enhanced. This occurs 
when an integer multiple of the energy scale of the 
tilting force Fd^ matches the energy difference between 
the initial state and the final state, i.e., approximately 
the band gap AE. This phenomenon has been ob- 
served experimentally :24. 25] and is called resonantly 
enhanced tunneling (RET). A deviation from the LZ 
prediction is expected in this case. In order to see 
whether the tunneling probability given by the standard 
LZ tunneling probability correctly predicts the height 
of a step corresponding to a single tunneling event, we 
fit a step function to our step-like survival probability 
and extract the height of each step. The result and the 
comparison to the LZ prediction is shown in the inset 




Fig. 2 The survival probability of the BEC in the ground 
band for a lattice depth Vo = 4.0 En, a constant Stark force 
Fo = 1.48, and various widths of the initial momentum dis- 
tribution of the BEC cloud, Ap/2p R = 0.05 (solid line); 
Ap/2p R = 0.2 (dashed line); Ap/2p R = 0.3 (dash-dotted 
line). Inset: step width At for several Ap (filled circles). The 
step width At is the distance between the two vertical dashed 
lines shown in the main panel for the case of Ap/2pB. = 0.05. 



of Figure [TJ When the system parameters are in the 
RET condition (e.g., I/7 ~ 1.05 and 1.6), the height 
of the steps of the survival probability increases and 
shows a significant deviation from the LZ prediction 
given in Eq. This behavior is seen in Figure [JJ 
where the survival probability for Fo = 1.48 decays 
much faster than for the other two cases. Since the 
derivation of Eq. (j4|) [27} does not take into account 
the actual Wannier-Stark level structure, which is 
necessary to describe the RET condition, Eq. (Q| or the 
survival probability derived from it cannot describe the 
enhancement of the tunneling probability due to RET. 

The other parameter which can affect the survival 
probability is the width Ap of the initial momentum 
distribution. By changing the trap frequencies to x and 
uj p one can prepare the initial distribution with different 
widths. Figure [2] demonstrates the dependence of the 
width of the steps At on Ap. The steps are smooth and 
partly washed out since the wave packet reaches the 
edge of the Brillouin zone earlier when it has a broader 
initial momentum distribution. Still the remnants of the 
steps cause a local deviation from the exponential decay 
of the tunneling probability. Nevertheless, the survival 
probability exhibits an exponential decay globally in 
time, i.e., on large time scales, for all the mentioned 
cases. 



2.1.2 Non-linear case - Gross-Pitaevskii equation In 
the regime of weak atom-atom interactions the effect of 
interactions is studied in the mean-field regime based on 
the Gross-Pitaevskii equation. The 3D Gross-Pitaevskii 
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equation can describe the dynamics of the entire Bose- 
Einstein condensate in terms of an equation of motion for 
a single particle wave function. Therefore, the following 
nonlinear equation describes the dynamics of interact- 
ing Bosc-Einstcin condensate atoms in a tilted periodic 
potential (V~2(x,t) = 0) with a = 1: 



ti 2 

H = + Krap(r) + Vi (x) + .93D |^(r,t) 2 | + Fx. 

(6) 

As mentioned in section [21 g^D is the coupling 
constant calculated from the s-wave scattering wave- 
length a s and the number of atoms in the conden- 
sate No. (^(r^) 2 ] is the local atomic density. As an 
estimate for the nonlinear term in Eq. (O, we define 

S3D|*'(r,t) 2 | 

C = — at the peak density of the initial 

state. The following effects can be seen by increasing 
the strength of a repulsive nonlinearity (<?3d > 0) in the 
system: 

— enhancement of the tunneling rate, 

— deviation from mono-exponential decay, 

— washed out steps (corresponding to damped Bloch 
oscillations). 

According to our results shown in Figure [3] for a RET 
case and the experimental results in [2~4"ll25j . the tem- 
poral behavior of the survival probability depends on 
the strength of atom-atom interactions. As can be found 
with more detail in [2~4"ll28| . the enhancement of decay 
rate is generic for repulsive interactions. The scaling of 
the decay rate as a function of nonlinearity is yet more 
interesting in the RET case (see the inset of Figure [3]). 
We can quantify the decay rate r of the survival prob- 
ability by globally fitting an exponential decay func- 
tion to the step-like curves of the survival probability. 
Such rates, for various nonlinearities, are depicted in the 
inset of Figure [21 A repulsive interaction initially en- 
hances the interband tunneling probability of the ultra- 
cold atoms |28j . Since the tunneling events occurring at 
different integer multiples of the Bloch period are corre- 
lated by the presence of the nonlinearity, a clear devia- 
tion from the mono-exponential decay is observed mak- 
ing the definition of a global decay rate r somewhat 
problematic [29]. A continuous change in the density of 
the condensate in time due to escaped particles from the 
system leads to a decreasing impact of the nonlinearity. 
Therefore, the time local rate of decay systematically 
decreases as the time increases. Additionally, the non- 
linearity leads to a dephasing and damping of the Bloch 
oscillations, not discussed here, but a discussion of this 
phenomenon can be found in (T.2 , 13 , 14 , 3T1I32J . 




i i i i i i i i i i i i i i i i i i i 

1 2 3 4 5 6 7 8 9 



Fig. 3 Time resolved survival probability in the ground 
band for the nonlinear case with Vb = 4.0 Er, Fo = 1.48 
(RET), a narrow width of the initial momentum distribution 
of the BEC cloud, Ap/2pn = 0.05, and composed of the 
following number of atoms prepared in a trap with u; x = 
2tt x 50 Hz and u p = 2tt x 100 Hz: iV = 2 x 10 4 , C » 0.2 
(open circles); No = 5 x 10 4 , C « 0.32 (filled squares); No = 
1 x 10 5 , C » 0.42 (open triangles) as compared to the linear 
case (?3d = 0, (solid line). Inset: Decay rate of the survival 
probability r at RET condition vs. No (filled circles). 



2.2 Impact of a time- dependent stochastic potential 
(noisy Wannier-Stark problem) 

Going back to the linear system, g 3 u = 0, the case of 
V2{x,t) ^ leads to a time-dependent stochastic poten- 
tial where we claim that we are able to control the dy- 
namics of the Bose-Einstein condensate atoms by chang- 
ing the characteristic parameters of the time-dependent 
stochastic phase <f>(t). a is no longer 1 and we can cal- 
culate it as defined above by Eq. (J3j, using the effective 
potential introduced below in Eq. ©. We use correlated 
noise for the time-dependent phase <p(t) for the second 
lattice. A standard example is exponentially correlated 
noise, which is characterized by a single correlation time. 
Such a noise can be obtained from linearly filtered white 
noise as follows 




where £ is a Gaussian white noise with zero mean and 
standard deviation equal to one. r is the correlation 
time of the noise and the strength of the noise is given 
by the parameter D. The introduced noise has a zero 
mean ((4>(t)) — 0) and an exponential correlation func- 
tion ((f>(t)(f)(s)) — D/rexp(— \t — s\ Jr). r describes the 
time-scale of the fluctuations and D/t its variance. The 
power spectrum of the exponentially correlated noise is 
given by a Lorentzian function as follows 

S (u) = jz - 5 „ . (8) 

V ^ 7r(l + 0J 2 T 2 ) ^ ' 
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Figure H] demonstrates the survival probability of the 
atoms in the ground band for various correlation times r 
and a fixed value of the variance D/t m 0.25 of the noise. 
The system parameters are Vq = 2.5 Er and Fq = 1.5 
(non-RET condition). The solid line shows the time evo- 
lution of the survival probability for the reference sys- 
tem introduced by Eq. (JTJ with g^ = 0, V2(x,t) = 
and a = 1. A nice step structure similar to the ones 
depicted in Figures [1] and [2] is observed for this case. 
The other curves in Figure U exhibit the time evolution 
of the survival probability for the temporally disordered 
system introduced in Eq. ((TJ) with — 0, V2(x,t) ^ 
and a ^ 1. As seen for all the cases the step structure 
is washed out. Considering the characteristic parameters 
of the exponentially correlated noise, different regimes of 
the noise are introduced. The regime of slowly varying 
noise is recovered when the noise has a large correla- 
tion time compared with the time scales of the system 
(e.g., Tb). The tunneling probability of the atoms in such 
a regime is suppressed (e.g., for r = 10 Tb shown by 
the dashed line in Figure U]) . Decreasing the correlation 
time the rate of tunneling increases but, nevertheless, 
for t > Tb the noise suppresses the tunneling compared 
with the reference system (see the data for r — 1 Tb 
as shown by the dot-dashed line in Figure 2]). For cor- 
relation times smaller than Tb, the noise recovers the 
regime of fast noise (e.g., for r = 0.05 Tb) and causes 
an enhancement in the tunneling rate (dot-dot-dashed 
line) compared with the reference system (thick solid 
line). Surprisingly a further decrease of the correlation 
time decreases the tunneling rate of the atoms (e.g., for 
t = 0.001 Tb depicted by the thin solid line). 

In order to understand the effect of the time- 
dependent stochastic potential on the system, we renor- 
malize the potential. The time-dependent potential can 
be replaced by a suitable static effective potential in the 
limit of small r. Such an effective potential can be cal- 
culated integrating over all possible phases giving: 



V e ft (x) = aV 



{TTx/d L ) + /3sm 2 Ux/dj) , (9) 



with a renormalization factor for the second lattice /3 



-2D It 



Eq. 



provides a time- independent potential 



which can be used to compute a as stated in Eq. (j3]). 
As seen in Figure U] the survival probability for a small 
r = 0.001 Tb (thin solid line) shows perfect agreement 
with the results achieved using the effective potential of 
Eq. © (filled circles). 

Keeping the noise parameters constant, we study the 
decay rate of the survival probability, T, of the conden- 
sate for various system parameters. The results of a scan 
over the Stark force are shown in Figure [Sfa). The simu- 
lations have been done for three values of the correlation 
time of the noise r = 0.0005,0.05,50 Tb, keeping the 
variance of the noise constant (here D/t 0.25). The 
three correlation times are chosen from the left shoulder, 
the peak point and the right shoulder of curves in Fig- 




3 4 5 6 7 8 9 10 



Fig. 4 Time evolution of the survival probability for Vo = 
2.5 E R , F Q = 1.5, and a « 0.82; The thick solid line is the 
decay of the reference system defined by Eq. JT|) with #30 = 0, 
V2(x,t) = and a = 1; in the presence of exponentially 
correlated noise with D/t m 0.25 and correlation times: r = 
10 T B (dashed line), r = 1 T B (dot-dashed line), r = 0.05 T B 
(dot-dot-dashed line), r = 0.001 Tb (thin solid line); results 
for the effective potential of Eq. (J9j) are shown by the filled 
circles. 



ure[5^6). The solid line depicts the decay rate of the refer- 
ence system (given by Eq.([rj with g 3 u — 0, V-2{x,t) = 
and a = 1). Noise in the system leads to a washing out of 
the RET peaks (present in the solid line of Figure [SJa)) 
in the decay rate, and, dependent on the correlation time 
of the noise, a suppression (for very large and very small 
values of r) or an enhancement of the tunneling (e.g., 
for r R3 0.05 Tb) can be obtained. 

We also ran a scan over the correlation time of the 
noise r, again keeping the variance of the noise constant 
(D/t ?s 0.25). The decay rate of the survival probability 
for three amounts of the Stark force (To = 0.95, 1.25, 1.5) 
is shown in Figure [5j6) . Comparing the symbols to the 
lines (which specify the decay rate in the corresponding 
reference system) one can realize that the decay rate is 
enhanced when the system parameters do not fulfill the 
RET condition, i.e., in this case for To = 0.95 and 1.5. 
For the case of To = 1.25 (RET) the tunneling rate is 
suppressed, and the symbols lie always below the ref- 
erence line (dot-dot-dashed line). The enhancement is 
pronounced in the range of t ph 0.005. ..0.2 Tb, corre- 
sponding to an energy scale equal or larger to/than the 
band gap AE. On the other hand, for small values of the 
correlation time, the noise recovers the regime of white 
noise, where the effective potential describes the dynam- 
ics of the system very well, and for large values of r, the 
system is in the regime of a slowly varying noise corre- 
sponding to energy scales that do not help to enhance 
the tunneling of atoms either. 
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3 Conclusion 

It is possible to control the Landau-Zener (LZ) tun- 
neling probability of the ultracold atoms from the 
ground band in tilted optical lattices. This control is 
possible by changing the system parameters such as the 
lattice depth and the Stark force, or by changing the 
initial condition which is given by the initial width of 
the momentum distribution of the BEC. All the men- 
tioned parameters can be easily tuned experimentally. 
Furthermore, our calculations showed that atom-atom 
interactions affect the LZ tunneling probability and a 
repulsive interaction typically leads to an enhancement 
of the LZ tunneling of the ultracold atoms from the 
ground band. Most interestingly, our results demon- 
strate that it is also possible to control the tunneling 
by adding noise to the system. By changing the noise 
parameters, the tunneling probability can be enhanced 
or suppressed. The noise can particularly enhance the 
tunneling probability when the system parameters are 
chosen far from the RET condition. 
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